The phase diagram in T - ft - N c space* 
Giorgio Torrieri , Stefano Lottini 6 , Igor Mishustin", Piero 

NlCOLINI a ' b 

a FIAS and b ITP, JW Goethe Universitat, Frankfurt, Germany 



We examine the phase diagram of hadronic matter when the number 
of colors, as well as temperature and density, are varied. We show that in 
this regime several phase transitions are possible, and we examine issues 
related to these transitions. 
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1. Introduction 

The large- N c limit [H [2] , while sharing several qualitative characteristics 
with QCD, exhibits significant simplifications. For this reason, it is a very 
popular technique for qualitative and semi-quantitative estimates of crucial 
QCD observables. For instance, the string-gauge correspondence is pro- 
foundly intertwined with the large N c limit, since one of the requirements 
for a calculable dual holographic description is large N c [3]. 

Yet, the existence of a smooth limit between N c = 3 and N c = oo is a 
conjecture which, for some observables, is known to be false. We know that, 
at zero chemical potential, N c = oo and finite Nf deconfinement has to be 
a phase transition. At N c = 3, Nf = 1,2 lattice QCD seems to indicate it is 
a cross-over [3]. At finite chemical potential, in the confining phase, nuclear 
matter is known to be a strongly bound crystal in the large N c limit [5]. It 
seems to be a liquid at iV c = 3. This means that, before predictions made at 
large N c can be made qualitatively reliable, the phase diagram in T — fi — N c 
space needs to be explored. The most current phenomenological application 
for these ideas is the existence of the so-called "quarkyonic matter" [6] . This 
talk explores these ideas, summarizing El [15] and drawing connections 
between these works. 
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2. The Van der Waals gas and large iV c 

Our first investigation of these issues [3] is to use the Van Der Waals 
(VdW) ansatz to see under what conditions can we obtain both the familiar 
nuclear matter phase diagram at N c = 3 and the commonly accepted large 
N c phase diagram at N c — > oo [6 1 . While the Van der Waals gas is certainly 
a very rough description, it is part of a universality class, to which the 
nuclear liquid-gas phase transition seems to belong. 

In the large N c limit, the only ^-invariant scale of the theory is Aqcd, 
the scale at which the 't Hooft coupling constant becomes A ~ O (1). While 
a precise value of this scale depends on the scheme used to calculate it, its 
roughly Aqcd ~ — 200 — 300 MeV [ID]. It is therefore natural to expect 
that any physical quantity is ~ f(N c )AQ CD , a dimensionless function of 
N c times a power of Aqcd set by the dimensionality d of the quantity. 
Henceforward we shall adopt this assumption, and, for brevity, set Aqcd 
to unity in the equations. In this notation, the Van Der Waals parameters 
a, b and the curvature correction become dimensionless a,/3,7 times the 
appropriate power of Aqcd (3 for a, 2 for /3,4 for 7), and the VdW equation 
becomes 

(p- x -<x)(p + Pp 2 - lf ?)=T (1) 

We start by noting that, naturally, a can only go to Aq CD at N c — Y 00, and 
is significantly larger than Aq CD at N c = 3. Therefore, 

(2) 

where Nn is a constant to be determined. The coefficients f3, 7 should, 
according to [2] [5] go as N c . Recent work [TTJ, however, has cast doubt on 
this assumption and proposed they go as ~ N® or ~ ln./V c . 

The chemical potential can be obtained by the textbook thermodynamic 
relation p = (dP/dp) T . Inverting, and writing in terms of \x q = /j,b/N c we 
have 

' IT f(p',T)dp' + F(T) 
Jo 



(3) 



where the first term is the nucleon mass and 



(4) 

p and P are the density at the phase transition, which could be liquid pi or 
gas p g (if the calculation is performed correctly the same chemical potential 
should come out). p^ g are in turn the solutions to the equation [JJ in the 
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region where this equation has two solutions (of which one is thermodynam- 
ically unstable). Obtaining all such solutions is trivial at the mathematical 
level through algebraically cumbersome. The reader can get the detailed 
results in [TJ. 

As can be seen (Fig. [TJ, to interpolate between the currenlty accepted 
nuclear liquid phase diagram at iV c = 3 (p c -C Aq CD ,T c -C Aqcd) and 
the "quarkyonic" phase diagram at N c — > oo (p c ~ Aq CD ,T c ~ Aqcd), 
Nn ~ O (10) must be used. How would we interpret physically such a 
value? Above all, what is the physical meaning of Nn? We conjecture 
that Nn tracks another dimensionless scale relevant at high density: The 
number of neighbors a nucleon has in a tightly packed nuclear material. 
AV, of course, is a function not of N c , but the (fixed) number of dimensions 
d and "packing scheme" . 

The sensitivity of the varying iV c equation of state to Njy is undestand- 
able due to the quantum nature of the problem and the fermionic nature of 
the quarks. The more neighbors, the more Pauli blocking of valence quarks 
must be important, and the more the presence of neighbors will disturb 
the configuration space part of the quark wavefunction inside the nucleons. 
Since, due to the uncertainty principle, any such disturbance of the nuclear 
wavefunction adds an energy of the order of the confinement scale ~ Aqcd, 
the nuclear repulsive core will be larger than the inverse of the nuclear 
separation up to the deconfinement temperature. If the number of colors 
is larger than Njy, this problem will not exist since it will be possible to 
arrange the color part of the wavefunction so the nearest quarks of neighbor- 
ing baryons will be of different colors. In this limit baryons can be tightly 
packed (interbaryonic separation ~ Aqcd) without the configuration space 
part of the baryonic wavefunction being disturbed. 

The number of neighbours scaling suggests that percolation phenomena 
might be relevant in the middle of the two regimes. In the next section we 
investigate this possibility. 

3. Percolation and dense matter at varying N c 

The key insight suggesting that interesting structures might be lurking 
in N c is that 3D bond-percolation exhibits a phase transition at compar- 
atively low critical link probability: for instance, p c ~ 0.25,0.18,0.12 for 
simple-cubic, body-centered-cubic and hexagonal-close-packed lattice, re- 
spectively [12]. Such values suggest that long-distance correlations on the 
quark level could occur even with a somewhat low percentage of quarks 
hopping between baryons, i.e. firmly in the confined phase. While below p c 
the characteristic correlation distance £ (~ cluster size) is ~ Aq CD , above 
the threshold this quantity explodes to the total system size in a comparable 
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Fig. 1. (color online) The phase diagram in T — fi space (left panel) and T — p space 
(right panel) as a function of N c (N c = 3,5,8,10,30,100, with increasing color 
corresponding to a line with higher T, /i, p) Top panels assume nuclear interactions 
~ N® or ~ lniV c , bottom panels as ~ N c 



amount of time. We leave the meaning of "correlation" vague, as it can be 
either a quark hop or a gluon exchange; in our context, it implies exchange 
of colour degrees of freedom within a confined tightly packed medium. We 
encode the likelihood of exchange between neighbouring baryons in a link 
probability p, to be compared with the percolation threshold p c in order to 
assess the formation of large-scale structures. 

Two baryons will be correlated if at least two quarks are correlated. 
One has to sum over all possible multi-quark configurations, resulting in a 
strong 7V c -dependence of p. We determine the latter by calculating the 
probability q = 1 — p that no exchanges happen between neighbouring 
baryons. Assuming the quarks inside the nucleon are uncorrelated (Fermi 
motion dominates), this probability factorises into a geometric distribution 
/a,b(x) for quarks to be at a certain (vector) position x, and a "squared 
propagator" transition amplitude F(d) for them: 

v = i - (q(d^ Nc)C ; (5) 

Q(i),ij = / /A(xi)dXi / f B (xj)dxj(l-F(\xi-Xj\)) . 



We assume a "hard-sphere" distribution for Ja,b (since we keep [i q fixed, 
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the distance between centres of neighbouring baryons is always 2A ( 



and a probability of exchange i«j based on a range of "reasonable" prop- 
agators, compatible with confinement (fast fall-out in configuration space 
at distances greater than ry ~ 1 in units of Kqcd) an d with the large N c 
limit of QCD, the interaction is ~ (? 2 ~ X/N c [HE]- The propagators we use 
are the simple 0-function in configuration space and the momentum-space 
0-function used in |13j . all normalised so their area is \vt/N c where rj- is 
the range of propagation (~ Aqq D ). In configuration space the transition 
amplitudes are, respectively, 



Other transition amplitudes, such as a Gaussian distribution in configura- 
tion space, were also tried with no significant modifications. This is un- 
surprising since, due to the fact that 3D percolation has a second-order 
phase transition at a certain p c < 1, the results we obtain below have some 
degree of universality: as long as the qualitative features of confinement 
are observed (the transition amplitude F{y) drops sharply above the scale 
vt, and the hadron density profile /a,b(x) has a central plateau of radius 
~ and a sharply decreasing tail outside), the results we show vary 

quantitatively but not qualitatively. 

The crucial parameter left is £ in Eq. One can easily see that £ = 1 
is in contradiction with the Skyrme crystal picture at large N c : in this 
picture, p(N c ) approaches a constant large- N c value from above: low N c 
nuclear matter would be more correlated (and hence more strongly bound) 
than high N c nuclear matter. Comparing strongly coupled N c — > oo nuclear 
matter [2] to the weakly bound nuclear liquid at JV C = 3 [ZJ, this is obviously 
not right. We therefore assume £ = 2 henceforward, natural if the link is 
actually realised by a gluon exchange rather than a quark flip (it is obvious 
this does indeed dominate at large N c ). 

In the large iV c limit for the case £ = 2, p asymptotically approaches 
unity. It is reasonable that this is the point where the "dense baryonic 
matter as a Skyrme crystal", theorised in [21 [6], is reached. If this is 
the case, however, one should remember that a percolation second-order 
phase transition occurs at a p c <C 1. Hence, keeping [i q ~ A-qcd fixed but 
varying iV c , the features of the Skyrme crystal should manifest not with 
a continuous approach, rather as a second order transition at a not too 
high N c , whose order parameter can be though to be the "giant cluster" 
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density. Below the critical N c there is little correlation between quarks of 
different baryons, while above this threshold they can correlate, with the 
distance boundary given only by causality. We reiterate that this is not 

1/2 

deconfinement since /i ~ h-QCD "C N c A-qcd independently of the number 
of colours, and the fraction of correlated quarks from different hadrons is still 
~ 0.1 — 0.3 <C 1 at the percolation transition. Right above this transition, 
therefore, the baryonic wavefunction should not be too different from the 
large- N c baryonic wavefunction described in [2]. The correlation distance of 
quarks will however be much larger than the baryon size. The features of 
this new phase are therefore similar to those of the quarkyonic matter [B] . 

Assuming the lowest 3D value p c = 0.12, appropriate for a closely packed 
hexagonal lattice, the critical number of colours is shown in Fig. [2] as a 
function of A and Tt- As can be seen, the critical number of colours is 
significantly larger than three for ~ ^qcdi ^ ~ 1- Considering Fig. [2] is 
a lower limit since p c is at its minimum (p c is significantly higher both in a 
Skyrme cubic crystal and in a disordered fluid), we can say that N c = 3 is 
disfavoured, although it can not be excluded. Changing temperature and 
\x q should further change the critical N c . Exploring this parameter space, 
and seeing how it relates to the confinement phase transition, is the subject 
of a forthcoming work. 

What are the phenomenological consequences of percolation? If by "cor- 
relation" we mean energy-momentum-exchange via quark tunneling between 
baryons, it is reasonable that pressure and entropy density ~ N c above the 
percolation threshold, while below it they stay ~ N®. This is because above 
the threshold, where interbaryon tunnelling is significant, "typical" excita- 
tions of the Fermi surface will be superpositions across baryons of baryon- 
localized quark- hole excitations, similar of conduction band electrons in a 
metal; while the localized excitation energy > Aqqd, the superposition 
makes its energy ~ Aqcd even if color degeneracy remains. Thus, the de- 
grees of freedom of the system above percolation will be delocalized weakly 
interacting quarks in a lattice of confining potentials, a picture compatible 
with [6]. Below the threshold, where tunnelling is negligible, excitations are 
either color singlets or of energy E S> Aqcd, suppressed below deconfme- 
ment. 



4. Confinement at finite N c and quantum gravity 

The transitions discussed here might also be visible with holographic 
techniques beyond the supergravity limit, since finite AA C _1 corresponds, 
in gauge/string duality, to the string coupling constant g s [3J. Thus, a 
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Fig. 2. (color online) Contour plot of the critical N c for the percolation transition in 
a hexagonally packed lattice as a function of the coupling A and range rx (in Aqcd 
units). The left panel assumes a 0- function correlation probability in configuration 
space, the right panel assumes a correlation probability based on the propagator 
used in [13]. Diagram covers 2 < N c < 80. See [8] for details 



transition at finite iV c will manifest itself as a transition between a classical 
supergravity regime and a quantum-gravitational regime. Since gravity to 
one loop still has not been solved for backgrounds dual to confining theories, 
we can not rigorously say anything further. 

However, in [9], we have made a first step in that direction using the 
reasonable Ansatz, developed earlier [HI US], that the effect of Quantum 
gravity corrections is to smoothen the Black hole singularity to a Gaussian 
of with £, the Plank length. The problem of the first order quantum black 
hole reduces, than, to the problem of the classical general relativity sourced 
by a Gaussian distribution of energy density p. The top energy-momentum 
tensor is then 

M ( r 2 \ 

with the other components of being fixed by isotropy, the classical 
Einstein equation 

1 3 

R/inu — ^Rdfiiy — 72 9{iv = SnGT^ (9) 

1 ADS 

and conservation laws ( Tf^ u = 0). Thermodynamics can be obtained by the 
usual relationship between the entropy s and the horizon area A, s = A/G, 
as well as Maxwell's relations |16j . 

We rely on [9] Q3] for details on solving this problem, and connecting it 
to the Hawking-Page phase transition [16] . and just give the result: 



F(r H ) = F hawking (r H ) + AF(r H ) 



(10) 
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where Fh aw king is the free-energy of a classical black hole in AdS space |16j 
and 



Once we examine the minimum structure of this free energy, we find a 
critical point: For a "critical" q = I/Lads — 0-18 there is a critical point 
after which the Hawking-Page phase transition becomes a cross-over (Fig. 
[3|). The resulting phase diagram (in T vs I) is shown to be in the universality 
class of the Van der Waals gas. 

The relevance of this paper to the topic at hand becomes apparent given 
that the Hawking-Page transition is thought to be dual to the deconfinement 
transition [3J. We know that at finite Nf/N c , deconfinement becomes a 
cross-over rather than a transition. Thus, this calculation can be thought 
of as an estimate of the critical Nf/N c for the deconfinement critical point 
in two-dimensions (since the gravitatational setup is in 3D). 

While the A0IS3 case presented here is considerably different from the 
AdS\ x S n or the AdS$ x S§ background examined in [3], our results fall 
neatly into the Van Der Waals universality class. Because of that, and the 
structure of the calculations in the previous section, we expect that the 
general structure of this phase diagram will be maintained in AdS\ x S n or 
any other background where a deconfinement phase transition exists. The 
critical q* will of course be different, but not the behavior around the critical 
point. This critical point will translate itself in the behaviour of N c of the 
boundary theory, once again in a way that is universal around the critical 
point but sensitive to the exact nature of the background theory. 

When the number of colors N c 3> the number of flavors Nj (the so-called 
T'Hooft limit), confinement is rigorously known to be a phase transition in 
all dimensions greater than one, since the symmetries of the system are 
expected to change. The low temperature phase will have acquired an extra 
Zfq symmetry, signifying a zero expectation value of the Polyakov loop 
[17j . In the high-temperature phase the Polyakov loop will aquire a finite 
expectation value, and Zfq will be spontaneusly broken. In the gravity 
picture, this can be thought of as a manifestation of the widely-believed 
cosmic censorship conjecture |18j . the idea that any singularity in spacetime 
is surrounded by an event horizon (and hence a discontinuity in the energy 
density at the boundary theory). Since the existance of the black hole 
tracks the confining Polyakov loop behaviour in [3J, making confinement 
into a cross-over would mean smoothening the black hole singularity. This 
can only be done, eg, by our ansatz, for q > 0. 

At finite number of flavors Nf, the Zfq symmetry is no longer exact at 
a fundamental level. Hence, it is natural to expect that at some critical 
jycmt " a cr iti ca l point" for deconfinement appears, where for higher N c 
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Fig. 3. The free energy F in units of the cosmological constant L for the system in 
equilibrium as a function of the Schwarzschild radius x+ from top to bottom for, 
in terms of q = i/L, q = 1/3 and q — 0.1 (solid lines). The dashed lines represent 
the critical case q — q* , while the dotted line is the classical Hawking-Page phase 
transition (q = 0). See [5] for details 



deconfinement is a phase transition and for lower N c it is a crossover. This 
is indeed what emerges from lattice simulations [3] in three dimensions with 
a finite number of flavors. The free energy with respect to the Polyakov 
loop expectation value (L) and energy E is therefore of the form 

F ((L) , E) = F„ ((L) ,E) + 6f(^-, (L) , . (11) 

Fgauge at the mixed phase is a function with two local minima, corresponding 
to the two phases. For high-enough Nf/N c , however, 5F could destroy one 
of the minima and make the phase transition into a crossover. The Van der 
Waals form of the free energy Eq. [TU] exactly parallels Eq. [TTl with the q 
parameter playing the role of N*/N c in a Gauge theory in 2+1 dimensions 
(since the bulk calculated here is 3+1 dimensions). 

5. Discussion and Conclusions 

The preceding three sections, at first sight, look like they deal with very 
different subjects. Indeed, the current status is that any correlation between 
them is conjectural. 

When N c is varied in the dense matter strongly coupled regime, devia- 
tions from the "classical" skyrme crystal arise due to the interplay between 
the number of colors and the number of neighbours iVjv in the densely packed 
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system. This "geometrical" picture is explored rigorously, in terms of an 
analytical model, in section [3] and shown to give a percolation-type phase 
transition. A similar behaviour can be argued in terms of the Pauli exclusion 
principle, and if it is assumed, as shown in section [21 leads to the transition 
between the loosely-bound nuclear liquid at N c = 3 and the tightly bound 
nucelar crystal at N c 3> Njy. Since the solid-liquid transition can be un- 
derstood as a percolation of crystal bonds [TI5], it is therefore compelling to 
identify the percolation transition with the physics described in section [51 
with the critical N c playing the role of "critical point" in T — \i — N c space 
(the transition in section [3] is necessarily second order). 

The connection with the Hawking-Page transition described in section [5] 
is, at the moment, not so clear. It is worth pointing out, however, that this 
transition is also in the Van der Waals universality class. While section|3]was 
at zero chemical potential (the black hole was not charged), the Hawking- 
page phase diagram for classical (£ = 0) charged black holes also falls into 
this class |20j. 

Thus, the phenomena examined at finite chemical potential in section 
[3] and [2] might be analogous to the dynamics examined in section H] at zero 
chemical potential. The closeness of the critical q* to the critical number 
of Nf/N c required for percolation in 2D [8] (6 colors per flavor) is in this 
sense intriguing. 

Physically, in AdS$, the Hawking-Page transition coincides with the 
transition from a gas of black holes to one large black hole. It is plausi- 
ble that the critical point to a Hawking-Page cross-over occurs because the 
black holes in the gas are actually linked by super-horizon quantum grav- 
itational interactions, modeled in our approach by "smearing" the super- 
horizon black hole distance by a Planck-sized width. In this scenario, the 
connection of percolation to the cross-over in deconfinement is straight- 
forward. Before this connection can be made more rigorous, however, a 
Gauge/string setup where percolation is relevant must be examined with a 
"smearing" analogous to section |H Work of this kind is in progress. 
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